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In this paper we study perturbatively an extension of the Stelle higher derivative gravity involving 
an infinite number of derivative terms. We know that the usual quadratic action is renormalizable 
but suffers of the unitarity problem because of the presence of a ghost (state of negative norm) in 
the theory. The new theory is instead ghost-free since the introduction of (in general) two entire 
functions in the model with the property do not introduce new poles in the propagator. The local 
high derivative theory is recovered expanding the entire functions to the lowest order in the mass 
scale of the theory. Any truncation of the entire functions gives rise to the unitarity violation but if 
we keep all the infinite series we do not fall into these troubles. The theory is renormalizable at one 
loop and finite from two loops on. Since only a finite number of graphs are divergent then the theory 
is super-renormalizable. We analyze the fractal properties of the theory at high energy showing a 
reduction of the spacetime dimension at short scales. Black hole spherical symmetric solutions are 
also studied omitting the high curvature corrections in the equation of motions. The solutions are 
regular and the classical singularity is replaced by a "de Sitter-like core" in r = 0. Black holes may 
show a "multi-horizon" structure depending on the value of the mass. 



I. INTRODUCTION TO THE THEORY 

One of the biggest problems in theoretical physics is 
to find a theory that is able to reconcile general relativ- 
ity and quantum mechanics. There are many reasons to 
believe that gravity has to be quantum, some of which 
are: the quantum nature of matter in the right-hand side 
of the Einstein equations, the singularities appearing in 
classical solutions of general relativity, etc. 

The action principle for gravity we are going to intro- 
duce in this paper is the result of a synthesis of minimal 
requirements: (i) classical solutions must be singularity- 
free, (ii) Einstein-Hilbert action should be a good approx- 
imation of the theory at an energy scale much smaller 
than the Planck mass, (iii) the spacetime dimension has 
to decrease with the energy in order to have a complete 
quantum gravity theory in the ultraviolet regime, (iv) the 
theory has to be perturbatively renormalizable at quan- 
tum level (this hypothesis is strongly related to the previ- 
ous one), (v) the theory has to be unitary, with no other 
degree of freedom than the graviton in the propagator, 
(vi) spacetime is a single continuum of space and time 
and in particular the Lorentz invariance is not broken. 
The last requirement is supported by recent observations. 

Now let us introduce the theory, step by step, starting 
from the perturbative non-renormalizable Einstein grav- 
ity, through high derivatives gravity theories (the Stelle 
theory of gravity will be our first example) onto the ac- 
tion which defines a complete quantum gravity theory. 
The impatient reader can skip to the end of the intro- 
duction for the candidate complete quantum gravity bare 
Lagrangian. 

Perturbative quantum gravity is the quantum theory 
of a spin two particle on a fixed (usually for simplicity 
is assumed to be flat) background. Starting from the 
Einstein-Hilbert Lagrangian 



(k^ = IGttGjv) we introduce a splitting of the inetric in 
a background part plus a fluctuation 
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then we expand the action in power of the graviton fluc- 
tuation /i'^" around the flxed background g'^'^ . Unlikely 
the quantum theory is divergent at two loops, producing 
a counter-term proportional to the Ricci tensor at the 
third power 
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In general, in d dimensions the superficial degree of diver- 
gence of a Feynman diagram is D = Ld + 2V — 21, where 
L is the number of loops, V is the number of vertices and 
/ the number of internal lines in the graph. Using the 
topological relation between V, I and L, L = I + 1 — V, 
we obtain D = 2 + {d — 2)L. In d = 4 the superficial de- 
gree of divergence D = 2 + 2L increases with the number 
of loops and thus we are forced to introduce an infinite 
number of higher derivative counter terms and then an 
infinite number of coupling constants, therefore making 
the theory not predictive. Schematically, we can relate 
the loop divergences in perturbative quantum gravity to 
the counter terms we have to introduce to regularize the 
theory. In short 
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n+2m=2+{d-2)L 



with "ri" and "tti" integer numbers. 

A first revolution in quantum gravity was introduced 
by Stelle [1| with the higher derivative theory 
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This theory is renormaUzablc but unfortunately contains 
a physical ghost (state of negative norm) , implying a vi- 
olation of unitarity in the theory: probability, described 
by the scattering S'-matrix, is no more conserved. The 
classical theory is unstable, since the dynamics can drive 
the system to become arbitrarily excited, and the Hamil- 
tonian constraint is unbounded from below. 

In this paper we generalize the Stelle theory to restore 
unitarity. This work is inspired by papers about a non- 
local extension of gauge theories introduced by Moffat, 
Cornish and their collaborators in the nineties Q]. The 
authors further extended the idea to gravity, having in 
mind the following logic [3|- They considered a modi- 
fication of the Feynman rules where the coupling con- 
stants {gi for electro-weak interactions and Gn for grav- 
ity) arc no longer constant but function of the momen- 
tum p. They checked the gauge invariance at all orders in 
gauge theory but only up to the second order in gravity. 
For particular choices oi gi{p) or Gn{p), the propagators 
do not show any other pole above the standard particle 
content of the theory, therefore the theory is unitary. On 
the other hand the theory is also finite if the coupling 
constants go sufficiently fast to zero in the ultra-violet 
limit. The problem with gravity is to find a covariant ac- 
tion that self-contains the properties mentioned before: 
finitencss and/or renormalizability and unitarity. 

In the rest of the paper we will mainly refer to the 
Tomboulis paper [J] that we have discovered only re- 
cently. The theory developed in [J| is very interesting, of 
great generality and mainly concentrated on gauge the- 
ories but in Q the author also concludes with an exten- 
sion of the idea to gravity. Here we will use the same 
notations as standard definitions. We arrive to the same 
conclusions in a different way, but the results are exactly 
the same. 

The action we are going to consider is a generalization 
of Stelle's theory 

5 = - Jv^{Rf,MOA)R'"' - RmA)R + JK-^R}. (6) 

where a{OA) and 13{D\) are now fixed functionals of the 
covariant D'Alembertian operator Da = D/A^ and A is 
a mass scale in the theory. 

In the remainder of this section we will summarize the 
steps and motivations that led us to the generalization 
^ and will conclude with an extended version ^. 

The covariant action (J6l) is a collection of terms, but 
the initial motivation when we started this project was 
the non local Barvinsky action Q 



^^C^" — i?pi/, 
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where Gf_,i, ~ R^i, — Rg^„/2 is the Einstein tensor. The 
action ([7]) looks like a non-local action, but indeed re- 
produces the Einstein gravity at the lower order in the 
curvature. Barvinsky has shown that the variation of ([7]), 
with respect to the metric, gives the following equations 



of motion (see later in this section for details) 



R,.. - ^9,.uR + 0{R%) - 0. 
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The next step is to modify ([7]) by introducing an extra 
operator in the action between the Einstein tensor and 
the Ricci tensor in order to have a well-defined theory at 
the quantum level without loss of covariance. 

We consider as guideline the value of the spectral di- 
mension in the definition of a finite and/or renormalizable 
theory of quantum gravity which, of course, is compatible 
with the Einstein gravity at low energy but also mani- 
fests a natural dimensional reduction at high energy. The 
dimensional reduction is of primary importance in order 
to have a quantum theory free of divergences in the ul- 
tra violet regime. The Stelle theory is characterized by a 
two dimensional behavior at high energy but this is not 
sufficient to have a well defined theory at the quantum 
level because of the presence of ghosts in the propagator. 
A first easy generalization of Q is the following action 
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where F{D/A'^) is a generic function of the covariant 
d'Alambcrtian operator which satisfies the already men- 
tioned properties that we are going to summarize below. 



(I) Classical limit 



lim F{a/A^) 
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if the limit is satisfied, the equations of motion are 
the Einstein equations plus corrections in i?^^ (we 
will explain this in more detail later in the paper). 

(II) Finiteness and/or renormalizability of the quan- 
tized theory. Our guiding principle is to find a well 
defined quantum theory and the dimensional reduc- 
tion of the spacetime at high energy. The Stelle 
theory, the Crane-Smolin theory [g, "asimptoti- 
cally safe guantur n gr avity" [3], "causal dynami- 
cal triangulation" [10|,"loop quantum gravity" [0] 
and "string theory" already manifest this prop erty 
with a high energy spectral dimension dg = 2 [l2|. 
However such reduction is insufficient if we want a 
unitary theory free from negative norm states. We 
can anticipate that, for the model that we are going 
to introduce in this paper, the spectral dimension 
is smaller than one in the ultraviolet regime. 

The general theory ^ was for the first time derived 
by Barvinsky [5| in the brane-world scenery and can be 
written in the following equivalent way 






d^X. 



-5<^ g Rf^iy- 



(11) 



In |5[ the author was interested in the infra-red modifica- 
tions to Einstein gravity. In the present paper, however. 



we are interested to ultra-violet modifications to grav- 
ity, as stressed in PH)) . The form pT|) of the action is 
particularly useful to highlight the classical limit and in 
the expression of the equations of motion. In analogy 
to the properties satisfied by the operator F(n/A^), the 
operator O has to satisfy the following limit in order to 
reproduce the classical theory (see the discussion related 
to formulas and ©) 



lim 0{D/A'^) = 1. 



(12) 



We now add some details about the Barvinsky derivation 
of the equation of motion. Taking the variation of (|lll) 
with respect to the metric we find (5| 
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Since SgR^^i, 
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+ V^e^ + V^e^ [3], integration 
by parts cancels the D operator at the denominator and 
the contribution of the gauge parameters e^ vanishes in 
view of the Bianchi identities, V^G^i^ = 0. All the com- 
mutators of covariant derivatives with the D operator 
in 0~^{Da)/0 give rise to curvature square operators. 
Also, the direct variation of the metric gives rise to cur- 
vature square terms. Then the equations of motion are 
very simple, if we omit the squared curvature terms: 
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If we truncate the theory to the linear part in the Ricci 
curvature, the equations of motions simplify to 
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where the second relation is a consequence of the Bianchi 
identities that wc impose on the solution. In particular, 
such relation implies that the conserved quantity has to 
be Sfii, = C'(n/A^)T^i^ in the truncation of the theory. 
These equations of motion are very interesting because 
black hole solutions, at least for a particular choice of the 
operator O, are not singular anymore, as recently showed 
in Q. 

At the classical level, equations p^ can be derived 
from the action ^ or (fTTj) but at the quantum level also 
the Einstein-Hilbert action and the high derivative terms 
introduced by Stelle are generated. Thus, the complete 
Lagrangian wc are going to study is 



A;R - p2 (Rn.R^" -Ir^]+ /?oi?' 
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-Rho (-Da) R - ^F^Lu{-al)F^ + C^M^.C 



R^, h2 (-Da) i?^" - \r h2 (-Da) r] (16) 



The operator D^ encapsulates the D'Alembertian of the 
fiat fixed background, whereas F^ is the gauge fixing 
function with the weight functional uj. In general we 
introduce two different functions /i2 and Hq. Those func- 
tions have not to be polynomial but entire functions with- 
out poles or essential singularities. While nonlocal ker- 
nels can lead to unitary problems, the functions ft,2 and ho 
do introduce an effective non-locality. However, since /12 
and ho are transcendental entire functions, their behav- 
ior is quite similar to polynomial functions and unitary 
problems do not occur. 

Let us assume for a moment that hi(x) = p„(x), where 
Pn{x) is a polynomial of degree n. In this case, as it will 
be evident in the next section, the propagator takes the 
following form 
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where we used the factorization theorem for polynomial 
and the partial fraction decomposition [J|. When multi- 
plying by k^ the left and right side of (|T7l) and considering 
the ultraviolet limit, we find that at least one of the co- 
efficients Ci is negative, therefore the theory contains a 
ghost in the spectrum. The conclusion is that /i2 and h^ 
can not be polynomial. 

We can also add to the action the Kretschmann scalar 
Rfj.iypaR'^'^'^'"^ but, for spacetimes topologically equivalent 
to the fiat space, we can use the Gauss-Bonnet topologi- 
cal invariant 

/ d^xV^[R^,p^R''''P^ - 4i?^,i?^'' + R^] (18) 

to rephrase the Kretschmann invariant in terms of R^ 
and R^ijR'^'^ already present in the action. 



II. BRST INVARIANCE AND GRAVITON 
PROPAGATOR 

We start by considering the quadratic expansion of the 
Lagrangian (|16|) in the graviton field fiuctuation without 
specifying the explicit form of the functionals /12 and /iq 
(if not necessary). Following the Stelle paper, we expand 
around the Minkowski background ry'"^ in power of the 
graviton field /i'^'^ defined in the following way 



^g^"" = if + nh^"". 



(19) 



The form of the propagator depends not only on the 
gauge choice but also on the definition of the gravita- 
tional fluctuation [8|. In the quantum theory the gauge 
choice is the familiar harmonic gauge d^h^^ = and 
the Green's functions are defined by the generating func- 
tional 






where Sg is the gravitational action defined in (J16p sub- alone. Other two BRST invariant quantities are 
tracted of the gauge and ghost terms and F"^ ~ F^^.h^'^ 



with F^^ 



SJ^di,. D'^ is the operator which gener- 



ates the gauge transformations in the graviton fiuctua- 
tion h^^^ . Given the infinitesimal coordinates transfor- 
mation x^^' ~ x^ + K^^, the graviton field transforms as 
follows 

+K{da,Ch"'' + d.^Ch"'" " Cdo.h^"' - d,,Chn- (21) 

We have also introduced a weighting gauge functional (to 
be precise the second term in (|20p '). which depends on a 
weight function uj{—0\) with the property to fall off at 
least like the entire functions /i2(fc^/A^), ho{k^/A'^) for 
large momenta [l| . 

When the gauge symmetry is broken by the addition of 
the gauge-fixing term, a residual transformation survives 
for the effective action which involves the gravitational, 
gauge-fixing and ghost actions terms. This is the BRST 
symmetry defined by the following transformation, which 
is appropriate for the gauge-fixing term, 

SbrstC'^ = -nC'^i-aDF^SX, (22) 

where 5X is a constant infinitesimal anticommuting pa- 
rameter. The first transformation in p2p is nothing but 
a gauge transformation generated by kC^SX, so the func- 
tional Sg is BRST invariant since it is a function of h^" 
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SBRSTidpC'C^) = 0, 
^ feRST(^r^") = 0. 



(23) 



The above transformation follows from the anticommut- 
ing nature of C" , i5A and the following commutation re- 
lation of two gauge transformations generated by ^'^ and 

77^, 






r^'^e)^KD^/{do.eri" - d^r,^ c)-m 



Given the second of (P^ . only the antighost C-r trans- 
forms under the BRST transformation to cancel the vari- 
ation of the gauge-fixing term. The entire effective action 
is BRST invariant 
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d4xa^;.i?rc^")=o. (25) 

Let us list the mass dimension of the fields in the gauge- 
fixed Lagrangian: [h^'^] = mass, [C^] = mass, [Cr] = 



mass, \K\ 



mass 



Now we Taylor-expand the gravitational part of the ac- 
tion (jlip to the second order in the gravitational pertur- 
bation h'^'^{x) to obtain the graviton propagator. In the 
momentum space, the action which is purely quadratic 
in the gravitational field, reads 



+{3 e [/3/2 - 3/3o«'fc' + iK'eh^ie/k^)] + 2r'^(fcVA')}p^",;,"Hfc) 

+k' [/3/2 - iPo^k' + 3K\'h,{k'/A')]{P^,-:\k) + V3[Pll--'^Hk) + P^"-r\k)]})h''^{k), (26) 



where we used the gauge F'^ = 9^/1^"^ and where we have 
introduced the projectors p^ 
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(27) 



where we defined the transverse and longitudinal pro- 



jectors for vector quantities 
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Using the orthogonality properties of ([23 
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we ean now 



invert the kinetic matrix in ((26)) and obtain the gravi- 
ton propagator. In the fohowing expression the graviton 
propagator is expressed in the momentum space accord- 
ing to the quadratic action 
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Let us consider the graviton propagator in the gauge ^ = 
0. In this particular gauge, only the first two terms in 
([211) survive. We will show in the next section that only 
the physical massless spin-2 pole occurs in the propagator 
when the theory is renormalized at a certain scale ^o- 
The renormalization group invariance preserves unitarity 
in the dressed physical propagator at any energy scale 
and no other physical pole emerges at any other scale. 



III. 



RENORMALIZABILITY 



In this section we want to find an upper bound to the 
divergences in quantum gravity; before doing this, we 
have to construct the entire functions /12 and ho- Look- 
ing at the first two gauge invariant terms in (|29p . we 
introduce the following notation 

h2{z) ^13- I32K^K^Z + K^K^Z h2iz), 

ho{z) ^ P- QPoK^k^z + Qk^K^z ho{z), (30) 

where z will be identified with — Da- 

Considering [J] , we require the following general prop- 
erties for the transcendental entire functions hi (i ~ 2,0): 

(i) hi(z) is real and positive on the real axis, it has 
no zeroes on the whole complex plane \z\ < -l-oo. 
This requirement implies that there are no gauge- 
invariant poles other than the transverse massless 
physical graviton pole. 

(ii) |/ii(z)| has the same asymptotic behavior along the 
real axis at ±00. 

(iii) There exists > such that 



hm Mz)\ ^ \zp 

|z|-i-l-oo 



7^2 
for the argument of z in the cones 



C = {z I - e < argz < -1-9 , tt - 9 < argz < n + 9}, 
for < 9 < 7r/2. 



This condition is necessary in order to achieve the 
supe-renormalizability of the theory. The necessary 
asymptotic behavior is imposed not only on the real 
axis, (ii) but also in conic regions surrounding the 
real axis. In an Euclidean spacetime, the condition 
(ii) is not strictly necessary if (iii) applies. 

Given the above properties, let us study the ultraviolet 
behavior of the quantum theory. From the property (iii) 
in the high energy regime, the propagator in momentum 
space goes as l/fc^'''+^ (see ([^H])). but also the n-graviton 
interaction has the same scaling, since it can be written 
in the following schematic way. 



-)■ /i" Drih {Or, + /i™ dhdy Drjh, 



(31) 



where h is the graviton field and hi is the entire function 
defined by the properties (i)-(ni). From ([^ . the superfi- 
cial degree of divergence (in four spacetime dimensions) 

is 

D = AL-{2-/ + 4)/ + (27 + A)V ^A- 2-/{L - 1). (32) 



In ([32)) we used again the topological relation between 
vertexcs V , internal lines / and number of loops L: I = 
V + L ~ 1. Thus, if 7 ^ 3, only 1-loop divergences exist 
and the theory is super-rcnormalizable. In this theory 
the quantities /3, 132, /3o and eventually the cosmological 
constant are renormalized, 

CKcn = C- V^l^i^^i? + \{Zx - 1) 

^p2{Z2 - iKRp.R'"' - ^R^) + MZo - l)R^}. (33) 

where all the coupling must be understood as renormal- 
ized at an energy scale fx. On the other hand, the func- 
tions hi are not renormalized. In order to better under- 
stand this point we can write the generic entire functions 
as series, hi{z) = J^r^'^rZ^- For r ^ 1 there are no 
countcrtcrms that renormahze a^ because of the superfi- 
cial degree of divergence ([321) • Only the coefficient qq is 
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rcnormalizcd but this is just a normalization convention. 
The non-trivial dependence of the entire functions hi on 
their argument is preserved at quantum level. 

Imposing the conditions (i)-(iii) we have the freedom 
to choose the following form for the functions hi, 



h2{z) 
ho{z) 



h2[z) — a + a2Z 
/lo(^) - a + a^z 



(34) 



for three general parameters a, a2 and ao- In general 
in order to be compatible with the conditions (i)-(iii), 
we have to take hi{z) — a cxp Hi(z), where Hi{z) is an 
entire function that exhibits logarithmic asymptotic be- 
havior in the conical region C. Since H{z) is an en- 
tire function. hi{z) has no zeros in all complex plane for 
\z\ < -|-oo. Furthermore, the non- locality in the action is 
actually a kind of non locality, because hi(z) is an expo- 
nential function and a Taylor expansion of hi{z) erases 
the denominator Da at any energy scale. For this reason, 
the entire functions hi(z) can be called quasi polynomial. 
The entire function H{z) which satisfies the property 
(iii), can be defined as 



fP-,+i{^) i_c(w) 



(35) 



where the following requirements have to be satisfied : 

a. p-y+i(z) is a real polynomial of degree 7 -I- 1 with 
P7+i(0)-0, 

b. ^(z) is an entire and real function on the real axis 
with C(0) = 1, 

c. \C{z)\ — >■ for \z\ — >■ 00 in the conical region C 
defined in (iii). 

Let us assume now that the theory is renormalized at 
some scale /ig. If we want that the bare propagator to 
possess no other gauge-invariant pole than the transverse 
physical graviton pole, we have to set 
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k2A2 
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(36) 



As pointed out in [Jl, the relations ([55)1 can be used to 
fix the introduced scale A in terms of the Planck scale 



If we fix 



/3(/xo) = a, 
PiilJ-o) 6/3o(Ato) 



(37) 



a2 ao 

the two mass scales are linked by the following relation 
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If the energy scale jjlq is taken as the renormalization 
point we find /12 = ^0 = f^i/J-o) expiJ(2;) := h{z) and 
only the physical massless spin-2 graviton pole occurs in 
the bare propagator. In the gauge ^ = the propagator 
in (|29p simplifies to 
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(27r)4 fc2 -H ie \^ a/i(fc2/A2) J 

- (2.)^a"(Pl^H-(^-)-4<-(^-)^ (^9) 

If we choose another renormalization scale fj,, then the 
bare propagator acquires poles; however, these poles can- 
cel in the dressed physical propagator because the shift 
in the bare part is cancelled with a corresponding shift 
in the self energy. This follows easily from the renor- 
malization group invariance. The same procedure is not 
applicable to the case ^2(2) = ^o(-z) = [l[, because the 
theory fails to be renormalizable when the unitarity re- 
quirement /32 = /?o = is imposed and an infinite tower 
of countcrtcrms has to be added to the action. 

An explicit example of h(z) = cxpiJ(z) that satisfies 
the properties (i)-(iii) can be easily constructed. There 
are of course many ways to choose C{z), but we focus 
here on the obvious exponential choice ({z) = exp(— 2;2)^ 
which satisfies property c. in a conical region C with 
G = it/A. The entire function H{z) is the result of the 
integral defined in 



Hiz) ^-[jE+r (0,p2^i(z))] + logb^+i(z)] , 



Re(p2_^i(z))>0, 



(40) 



where je = 0.577216 is the Euler's constant and 
r(a, z) = J^ t°~^e~*dt is the incomplete gamma func- 
tion. If we choose p-y+i(z) = z^~^-^, H{z) simplifies to: 



Rc(z2''+2) > 0. 



(41) 



Another equivalent expression for the entire function 
H{z) is given by the following series 



+00 



n=l 

Re(p2^i(z)) > 0. 
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(42) 



For p^+i(z) = z'^^^ the O angle, which defines the cone 
C, is 8 = 7r/(47 4-4). According to the above expression 
(|42|) we find the following behavior near z = Q for the 
particular choice p-y+i(z) = 2'*'+^, 



H{z) 



^27+2 



Z 



47+4 
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267+6 
36 



0(z6^+^), (43) 



where the Taylor expansion is the exact one for Arg(2;) < 
7r/8 but wc already have a stronger constraint on the 




FIG. 1: Plot of \h{z)f for z real and a = 02 = 1 (solid line). 
In the first plot z e [—10, 10] and in the second one z £ [—2, 2]. 
The dashed line represents the asymptotic limit for large real 
positive and negative values of z. The asymptotic behavior is 
\h{z)f ^1.8 z'^. 



cone 
lim._ 



C, Arg(z) < 



Q ~ 7r/(47 + 4). In particular 



IV. SPECTRAL DIMENSION 

In this section, we calculate the spectral dimension of 
the spacetime at short distances, showing that the renor- 
nializability, together with the unitarity of the theory, 
implies a spectral dimension smaller than one. Let us 
summarize the definition of spectral dimension in quan- 
tum gravity. The definition of spectral dimension is bor- 
rowed from the theory of diffusion processes on fractals 
|ll| and easily adapted to the quantum gravity context. 
Let us study the Brownian motion of a test particle mov- 
ing on a d-dimensional Riemannian manifold M with a 
fixed smooth metric g^u{x). The probability density for 
the particle to diffuse from x' to x during the fictitious 
time T is the heat- kernel Kg{x, x';T), which satisfies the 



^ We can do a simple choice of the entire function H{z), which 
gives rise to a condition stronger than (iii). If wc take H{z) = 2^, 



then 



h.2{z) = ho{z) 



(44) 



where again 2 = — D/A^. Another possible choice we wish to 
analyze is h2 = ho = ae^, because of its connection with the 
regular Nicolini-Spallucci black holes [lil . 

For the functions given in l|44|l . the upper bound in l|32|l can be 
derived as a truncation of the exponentials and the result does 
not change: we have divergences at one loop, but the theory is fi- 
nite for L > 1. The exponentials in 1441 1 improve the convergence 
properties of the theory and the propagator is: 



-tv^fypfT yfi) ■ 



2e 



-fe^A-* 



(27r)'l a(fc2-|-ie) 



,(2) 



ik)-2PlSp^ik)]. (45) 



For hi{z) = aexp(2) the exponential in the above propagator is 
replaced with exp(— fc-^/A^). 



heat equation 



dTKg{x,x';T) 



AfKg{x,x';T) 



(46) 



where A"'^ denotes the effective covariant Laplacian. It is 
the usual covariant Laplacian at low energy but it can un- 
dergo strong modification in the ultra-violet regime. In 
particular, we will be interested in the effective Lapla- 
cian at high energy and in relation to the flat back- 
ground. The heat-kernel is a matrix element of the 
operator exp(TAg), acting on the real Hilbert space 
L'^i-M, ^/gd'^x), between position cigenstatcs 



Kg{x,x';T)^{x'\cxp{TAf)\x). 



(47) 



Its trace per unit volume. 



Pg{T) = y-1 fd^x^^Kg{x,x;T) 

= V-'^TTexp{TAf) (48) 

has the interpretation of an average return probability. 
Here V = J df^x ^ denotes the total volume. It is well 
known that Pg{T) possesses an asymptotic expansion for 
T ^ of the form Pg{T) = {AttT)-'^''^ Y.n=o ^" ^"- The 
coefficients A^ have a geometric meaning, i.e. Aq is the 
volume of the manifold and, if d = 2, Ai is then propor- 
tional to the Euler characteristic. From the knowledge of 
the function Pg (T) one can recover the dimensionality of 
the manifold as the limit for small T of 



^ _,dlnP,(T) 
dlnT 



(49) 



If we consider arbitrary fictitious times T, this quantity 
might depend on the scale we are probing. Formula (j49p 
is the definition of fractal dimension we will use. 

From the bare graviton propagator (j39p we can eas- 
ily obtain the heat-kernel and then the spectral dimen- 
sion of the quantum spacetime. In short, in the momen- 
tum space the graviton propagator, omitting the tenso- 
rial structure that does not affect the spectral dimension, 
reads 



D{k) oc 



1 



F/l(fc2/A2)' 



(50) 



We also know that the propagator (in the coordinate 
space) and the heat-kernel are related by |la | 



G{x,x') 



+ C» 



dTKg{x,x';T) 



/r+00 
d^fce*'^^^-^') / dTKgik;T), (51) 



where G'(x,a;') cx/d^fc exp[ik{x — x')]D{k) is the Fourier 
transform of (|50|) . Given the propagator ([50|) . it is easy to 
invert (|5ip with the heat-kernel in the momentum space. 



Kg{k-T) (xcxp[-fc2/i(fcVA^)r], 



(52) 



which is the solution of the heat-kernel equation p6)) with 
the effective operator 



Af = M-A,/A2)A,, 



(53) 



which goes like {—Ag)''^'^ Ag at high energy. The neces- 
sary trace to calculate the average return probability is 
obtained from the Fourier transform of ([5^ . 



Kg{x,x';T)Gc /d4fce-'='''(fcVA=)Tg^fc(--x'). (54) 



Now we are ready to calculate the average return proba- 
bility defined in 



(55) 



Pg{T)o, d'ke-^^-''^'^^^'^^. 




FIG. 2: Plot of the spectral dimension as a function of the 
fictitious time T for the special case 7 = 3 in (|42p . The 
lowest value in the picture is ds — 4/5 at high energy, but at 
low energy the spectral dimension flows to ds — 4. 



From the requirement (iii) we know that, at high energy, 
h{k) ^ k'^'' and then h{k) ^ fc^'*'+^; therefore, we can cal- 
culate the integral ([55]) and then the spectral dimension 
defined in (1491) for small T will be 



(56) 



PgiT)^ 



T2+^ 



7 + 2 



The parameter 7^8 implies that the spectral dimension 
is ds < 1, manifesting a fractal nature of the spacetime 
at high energy. 

We can calculate the spectral dimension at all energy 
scales as a function of the fictitious time T using the 
explicit form of the entire function H{k^ / h?) given in 
([42|) . Integrating numerically ([55|). we can plot directly 
the spectral dimension achieving the graphical result in 

FigHF. 



^ For the operators introduced in the previous section 
exp(— n/A^)" (n = 1,2), the propagator scales as 



D{k) oc 



fc2 



(57) 



and the spectral dimension goes to zero at high energy. In par- 
ticular, for n = 1 the heat-kernel can be calculated analytically. 



K{x,x';T) 



e~4CT + l/A2) 

[47r(T-Hl/A2)]2 



(58) 



as it is easy to verify by going back to the propagator JSTJ . Now, 
employing Eq. ll49l l. we find that the spectral dimension is 



4T 



1/A2 



(59) 



which clearly goes to zero for T — > and approaches ds = 4 for 

T -s> -1-00. 



V. BLACK HOLES 

In this section, we want to solve the equation of mo- 
tion coming from the renormalized theory in the case of 
a spherically symmetric spacetime. Let us start from the 
classical Lagrangian that we rewrite rearranging the pa- 
rameters in a different way, 



C 






(/?0 






P2 



a2 



1 



6k2A2 



R^ -a C 



i?„,.},(60) 



>(-nA)-i 



K2n 



i-liv I 



where h{z) := expi7(z)^. From [5[, the equations of 
inotion for the above theory up to square curvature terms 
are 

G^, + 0{Rl,) + 0(V^V,i?p,) = SttGat/i-IT^,, (61) 

where we omitted the argument of /i(— Da). Since we are 
going to solve the Einstein equations neglecting curvature 
square terms, then we have to impose the conservation 
V'^(/i^^r^,y) = in order for the theory to be compati- 
ble with the Bianchi identities. For the exact equations 
of motion the Bianchi identities are of course satisfied 



^ In general a differential equation with an infinity number of 
derivative has not a well-defined initial value problem and it 
needs an infinite number of initial conditions. It is shown in 
Il8l that in a general framework each pole of the propagator 
contributes two initial data to the final solution. This is pre- 
cisely our case because the only pole in the bare propagator is 
the massless graviton and the theory has a well defined Cauchy 
problem. 



because of the diffeomorphisms invariance. The condi- 
tion V^{h~^T^^) = compensates the truncation in the 
modified Einstein equations (|6T|) . 

Our main purpose is to solve the field equations by 
assuming a static source, which means that the four- 
velocity field u^ has only a non- vanishing time-like com- 
ponent u^ = (^",0) w" = (— ff^")"^/^. We consider 
the component T° o of the energy-momentum tensor for 
a static source of mass m in polar coordinates to be 
TOq = -m(5(r)/47rr2 [lll,[l3l'^- The metric of our space- 
time is assumed to be given by the usual static, spheri- 
cally symmetric Schwarzschild form 



ds^ = -F{r)di 



2 , °^^ , ^2^2 



F{r) 



F{r) = 1 



2m{r) 



(62) 



The effective energy density and pressures are defined by 



h-'T'',. 



SttG 



N 



Diag(-p^P;,PI,P^ 



(63) 



For later convenience, we temporarily adopt free-falling 
Cartesian- like coordinates [13|, we calculate the effective 
energy density assuming p^+i(z) = z^ in H{z), 

d^k 



^p sin(p) dp, (64) 



(27r)3 
2m '■+°° 



(2^) 






where r = |x| is the radial coordinate. Here we intro- 
duced the Fourier transform for the Dirac delta function 
and we we also introduced a new dimensionless variable 
in the momentum space, p = kr, where k is the physical 
momentum. The energy density distribution defined in 
(|64|) respects spherical symmetry. We evaluated numeri- 
cally the integral in (|64|) and the resulting energy density 
is plotted in FigUl In the low energy limit we can ex- 
pand H{z) for z = — D/A^ <C 1 and we can integrate 



^ Usually, in General Relativity textbooks, the Schwarzschild so- 
lution is introduced without mentioning the presence of a point- 
like source. Once the Einstein equations are solved in the vac- 
uum, the integration constant is determined by matching the so- 
lution with the Newtonian field outside a spherically symmetric 
mass distribution. Definitely, this is not the most straightfor- 
ward way to expose students to one of the most fundamental 
solutions of the Einstein equations. Moreover, the presence of a 
curvature singularity in the origin, where from the very begin- 
ning a nite mass-energy is squeezed into a zero- volume point, is 
introduced as a shocking, unexpected result. Against this back- 
ground, we show that, once quantum delocalization of the source 
is accounted, all these flaws disappear. From this follows that 
for us there is only one physical vacuum solution and this is the 
Minkowski metric. In other words, the Schwarzschild metric is a 
vacuum solution with the free integration m equal to zero. 



analytically (|64|) 



(2^) 






+00 



'p sin(p) dp. (65) 



The result is really involved and the plot is given in FiglSj 
however, the Taylor expansion near r w gives a con- 
stant leading order 



p%r) 



m 



A3 



32 2Vi6r(ll/l6)r(7/8)r(5/4) 



+ 0[r'^). (66) 




FIG. 3: Plot of the energy density for m = 10 in Planck units 
assuming A = mp. The solid line is a plot of (|64|) without any 
approximation, the dashed line refers to the energy density 
profile ([651) in the limit -D/A^ < 1. 



The covariant conservation and the additional condi- 
tion, goo = —grr: Completely specify the form of h^'^T'^^ 
and the Einstein's equations reads 

dm{r) 2 



dr 
1 dF 
F~d^ 



«^3^ 



2(m(r) +ATTP^r 
r{r — 2m{r)) 



dr ^ 2F dr ^^ 



-{Pl-Prl, (67) 

r 



Because of the complicated energy density profile, it is 
not easy to integrate the first Einstein equation in (|67p 



m{r) =4:71 dr'r''^ P^ir'). 
Jo 



(68) 



However, the energy density goes to zero at infinity, re- 
producing the asymptotic Schwarzschild spacetime with 
m{r) w m (constant). On the other hand, it is easy to 
calculate the energy density profile close to r « since 
H{z) -^ logz'^ for z -> -1-00 (or r -> in ([Ml))- In this 
regime mir) ex m A^r^ and for a more general monomial 
p-y+iiz) = zT+\ m{r) ex m{Ar)'^''+^ . The function F{r) 
in the metric, close to r w 0, is 



F(r) Si 1 — cm 



^27+2^27+1 



(69) 



where c is a dimensionless constant. 
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We show now that the metric has at least two horizons, 
an event horizon and a Cauchy horizon. The metric inter- 
polates two asymptotic flat regions, one at infinity and 
the other in r = 0, so that we can write the g~j^ = F 
component in the following way 



Fir) 



1 



2m/(r) 



(70) 



where /(r) — > 1 for r — > oo, /(r) ex r^''^^ for r — > and 
/(r) does not depend on the mass m. The function F{r) 
goes to "1" in both limits (for r -^ +oo and r — > 0) and, 
since tti is a multiplicative constant, we can always choose 
the mass m for a fixed value of the radial coordinate r, 
such that F{r) becomes negative everywhere. From this 
it follows that the function F{r) must change sign at least 
twice. The second equation in (|67|) is solved by P^ = —p'^ 
and the third one defines the transversal pressure once 
known the energy density p"^. Given the lapse function 
F{r) in (|69l) . we can calculate the Ricci scalar and the 
Kretschmann invariant 



R = cm A2T+2 (27 + 2)(27 + 3) r^''-^ 

= 4 c^ m^ A4t+4 (47^ + 47=^ + 57^ + 47 + 2) ' 



(71) 



.47 



By evaluating the above curvature tensors at the origin 
one finds that they are finite for 7 > 1/2 and in particular 
for the minimal super-renormalizable theory with 7^3. 

The entire function h{z) is able to tame the curvature 
singularity of the Schwarzschild solution at least for the 
truncation of the theory here analyzed. We think that 
the higher order corrections to the Einstein equation will 
not change the remarkable feature of the solutions found 
in this section. 

Besides the analysis exposed above, we can integrate 
numerically the modified Einstein equation of motions 
(j6ip for the two energy densities defined respectively in 
((64|) and (j65|) . Using the integral form of the mass func- 
tion (|68p . wc achieve the metric component F{r) defined 
in ((62|) . The numerical results are plotted in FigH] and 
FigOfor different values of the ADM mass m. The met- 
ric function F(r) can intersect zero times, twice or more 
than twice the horizontal axis relative to the value of the 
ADM mass m. This opens the possibility to have ^^multi- 
horizon black holes^^ as an exact solution of the equation 
of motions (jlSp . 



VI. STRUCTURE OF THE INTERACTIONS 

We have already shown that the theory is well-defined 
and power-counting super-renormalizable. However, the 
calculations are not easy beyond the second order in the 
graviton expansion. In this section, we give a sketch of 
how to proceed in the graviton expansion. The reason 
for the plots in FiglHlis to give an operative definition of 
exp H{z). One possible approximation in the interactions 




20 40 60 80 100 




FIG. 4: The two plots show the function F(r), assuming the 
infrared energy density profile ((64]). The two plots differ in 
the value of the ADM mass, which is tti = 1 in the fist plot 
and 771 = 10 in the second plot (in Planck units, assuming 
the fundamental scale A to be the Planck mass). A crucial 
property of those black holes is the possibility to have "multi- 
horizon black holes" depending on the mass value. For m — 
10, for example, we have six horizons according to the second 
plot. 



is the following replacement in the graviton expansion for 
the minimal rcnormalizable theory with 7 = 3, 



.H{z} 




72 



288 



for z > 1.3 , 
for z< 1.3. 



(72) 



At tree level, or loop amplitudes at high energy, we can 
just replace expH{z) with z^, at low energy with the 
second expansion defined in ([72)) for z < 1.3 and proceed 
in the calculation, gluing together the results in the two 
different regimes. 

Let us recall again the classical Lagrangian, 



C 



P2 



a2 



k2A2 



1 



{RpiyR^" — t;R ) 



I3q 



ao 



6k2A2 



i?2. 



aG^"" 






■Rf^i^j, 



(73) 
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20 40 60 80 100 





FIG. 6: The first plot is the function exp H(z) for p-y+i = P4 = 
2*. The solid line represents the exact function, the large and 
small dashed lines represent the same function for large and 
small value of z respectively: z^ and exp(z^''/36 — z^^ /8 + 
z^/2). In the second plot, the small dashed line represents 
a further simplification of the same function: expH{z) « 



1 + 272 ■ 



772 + 2-^7288. 



tree level, n-points functions will be obtain by an inter- 
polation of the amplitude calculated in the two different 
regimes A: < A and fc > A, using respectively the £J^ 
and >Cuv- ^^ 1*-"°? amplitudes, we should integrate the 
interactions terms coming from C^^ up to fc < A and the 
interactions coming from ^Cyy in the range A < fc < +00 
in the same amplitude. 



FIG. 5: The first plot shows the function F{r) for the energy 
profile (|65|) and H{z) defined in (|4H) with the parameter 7 — 
3. The ADM mass values are m — 1 and tti = 10 (in Planck 
units) for the first and the second plot respectively. 



where exp H{z) is defined in (|72[) . This Lagrangian inter- 
polates between the Einstein-Hilbert Lagrangian at low 
energy and a high energy theory living in a spacetimc of 
spectral dimension ds ~ 4/(7 -h 2). 

A first approximation (but also an operative way to 
proceed) is to replace the interaction Lagrangian in the 
UV with the following truncation 



/'int 



a 



UV 



2^2 



^gC^" 



-a 



A2 



Rf.^ for fc > A, (74) 



(fc is the energy scale) and the infra-red Lagrangian with 



/•int 



+ 



IR 



{l3o 



^) (i^^.i?^'^ - ^R') 



ao 



6«;2A2 






for fc < A. In (fM]) and ([75]) we used the same expansion 
of ((72|) but for a general value of 7. On the other hand, 
the propagator is the same in both regimes and is given 
in (j29p . The three graviton interactions can be obtained 
performing an h'^'^ power expansion in (|74p and ()75p . At 



VII. CONCLUSIONS 

In this paper we studied a new perturbative quantum 
gravity theory with a "gentle non-local character." We 
have shown that it is possible to build a perfectly well- 
defined quantum gravity theory without extra poles in 
the graviton propagator. 

The properties required for the theory introduced and 
studied in this paper were the following: 

(i) the theory should reproduce general relativity in 
the infra-red limit; 

(ii) black hole solutions of the classical theory have to 
be singularity free; 

(iii) the theory should be perturbatively renormalizable 
or super-renormalizable or finite; 

(iv) the spectral dimension should decrease at short dis- 
tances; 

(v) the theory has to be unitary with no other degrees 
of freedom than the graviton. 

All the above properties are satisfied by our quite re- 
strictive class of actions, which differ uniquely for the 
choice of an entire function H{z). This class of theo- 
ries can not be renormalizable or finite but such theo- 
ries turn out to be super-renormalizable, since only one 
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loop Fcynman diagrams diverge, implying a renormal- 
ization of just three coupling constants. The propagator 
has only one pole in the graviton mass shell. The min- 
imal super-renormalizable theory we built has spectral 
dimension ds = 4/5 in the ultra-violet regime and is four 
dimensional in the infrared limit. 

We have also considered a truncation of the classical 
theory, showing that spherically symmetric black hole so- 
lutions are singularity- free. As for the solutions in |19j . 
we have black holes only if the mass is bigger than the 
Planck mass (if we assume the fundamental scale in the 
theory to be the Planck mass). The new black hole so- 
lutions are more properly ^^multi-horizon black holes,^ 
showing a very reach spacetime structure depending on 
the value of the ADM mass. 

Future work will focus on the following subjects: 

• the cosmological singularity problem at the classical or 
semiclassical level. Some preliminary work has been 
already done in [201 ! 

• the high curvature corrections to the black hole solu- 
tions presented in this paper. 

Future work will also take the following directions: 

o we will consider more in detail the connection between 
non-locality and the fractality of the spacetime (see in 



particular the recent work by Calcagni in |T2|): 

o we will reconsider super-renormalizable gauge theories, 
with particular attention to the grand-unification of 
the fundamental interactions with or without gravity. 
It this work, we will try to put together ideas coming 
from [4| and the more recent paper |21| . 

Possible simplifications of the theory are: 

* a scalar field theory with the same non-local structure; 

* a simplification of the metric to the conformal form 
g^ij = QSx) rj^i, and therefore a quantization of the con- 
formal factor Vtix). 
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